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Abstract 

We study a model of random surfaces arising in the dimer model on the honeycomb 
lattice. For a fixed "wire frame" boundary condition, as the lattice spacing e —* 0, 
Cohn, Kenyon and Propp [3J showed the almost sure convergence of a random surface 
to a non-random limit shape So- In [11], Okounkov and the author showed how to 
parametrize the limit shapes in terms of analytic functions, in particular constructing 
a natural conformal structure on them. We show here that when So has no facets, 
for a family of boundary conditions approximating the wire frame, the large-scale 
surface fluctuations (height fluctuations) about So converge as e — ► to a Gaussian 
free field for the above conformal structure. We also show that the local statistics of 
the fluctuations near a given point x are, as conjectured in [3], given by the unique 
ergodic Gibbs measure (on plane configurations) whose slope is the slope of the tangent 
plane of So at x. 
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1 Introduction 

1.1 Dimers and surfaces 

A dimer covering, or perfect matching, of a finite graph is a set of edges covering all 
the vertices exactly once. The dimer model is the study of random dimer coverings 
of a graph. Here we shall for the most part deal with the uniform measure on dimer 
coverings. 

In this paper we study the dimer model on the honeycomb lattice (the periodic 
planar graph whose faces are regular hexagons), or rather, on large pieces of it. This 
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Figure 1: Honeycomb dimers (solid) and the corresponding "lozenge" tiling (green). 



model, and more generally dimer models on other periodic bipartite planar graphs, are 
statistical mechanical models for discrete random interfaces. Part of their interest lies 
in the conformal invariance properties of their scaling limits [HJ [9] . 

Dimer coverings of the honeycomb graph are dual to tilings with 60° rhombi, also 
known as lozenges, see Figure [TJ Lozenge tilings can in turn be viewed as orthogonal 
projections onto the plane P\\\ = {x + y + z = 0} of stepped surfaces which are 
polygonal surfaces in M 3 whose faces are squares in the 2-skeleton of Z 3 (the stepped 
surfaces are monotone in the sense that the projection is injective), see Figures ll|31 
Each stepped surface is the graph of a function, the normalized height function, 
on the underlying tiling, which is linear on each tile. This function is defined simply 
as v3 times the distance from the surface to the plane Pm. (The scaling factor y/3 is 
just to make the function integer- valued on Z 3 .) 

1.2 Results 

We are interested in studying the scaling limit of the honeycomb dimer model, that 
is, the limiting behavior of a uniform random dimer covering of a fixed plane region U 
when the lattice spacing e goes to zero. Equivalently, we take stepped surfaces in eZ 3 
and let e — > 0. As boundary conditions we are interested in stepped surfaces spanning a 
"wire frame" which is a simple closed polygonal path 7 e in eZ 3 . We take j e converging 
as e — * to a smooth path 7 which projects to dU. 

1.2.1 Limit shape 

Let U be a domain in Pm, and 7 be a smooth closed curve in M 3 , projecting orthog- 
onally to dU. For each e > let 7 e be a nearest-neighbor path in eZ 3 approximating 
7 (in the Hausdorff metric) and which can be spanned by a monotone stepped surface 
S e , monotone in the sense that it projects injectively to Pm, or in other words it is the 
graph of a function on Pm. See for example Figure [3] (although there the boundary is 
only piecewise smooth). 
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The existence of such an approximating sequence imposes constraints on 7, see 
[3j [5], as follows: the curve 7 can be spanned by a surface E, which is the graph of 
a continuous function on U, and such that the normal v to the surface points into 
the positive orthant M> - Conversely, any such curve 7 can be approximated by 7 e as 
above, see [3]. The condition of positivity of the normal to S can be stated in terms 
of the gradient of the function h whose graph is X: this gradient must lie in a certain 
triangle. The formulation in terms of the normal is more symmetric, however. 

For a given j e there are, typically, many spanning surfaces T, e and we study the 
limiting properties of the uniform measure on the set of T, t as e — > 0. 

For a surface T, t spanning j e , let h t : P%\\ — > R be the normalized height function, 
defined on the region enclosed by U e := 7rm(7 e ), whose graph is S e . 

Under the above hypotheses Cohn, Kenyon and Propp proved the existence of a 
limit shape: 

Theorem 1.1 (Cohn, Kenyon, Propp |3j) The distribution ofh e converges as e — > 
a.s. to a nonrandom function h : U — > R. T/ie function h is the unique function h 
which minimizes the "surface tension" functional 

min / a(Vh) dx dy, 
h Ju 

where, in terms of the normal vector (p a ,Pb,Pc) S R 3 to the graph of h scaled so 
that p a + Pb + Pc = 1, we /icwe a(p a ,p b ,p c ) = -±(L(irp a ) + L(7rp 6 ) + L(irp c )) and 
L(x) = — Jq 1 log 2 sin t is the Lobachevsky function. 

Here the minimum is over Lipschitz functions whose graph has normal with non- 
negative coordinates. Equivalently, these are functions whose gradient lies in a certain 
triangle. In the above formula the surface tension a is the negative of the exponential 
growth rate of the number of discrete surfaces of average slope Vh. 

The function h is called the asymptotic height function. Its graph is a surface 
£0 spanning 7. 

1.2.2 Local statistics 

Suppose that the gradient of h is not maximal at any point in U, i.e. the normal to So 
has nonzero coordinates at every point of U. In this paper we show that, if the precise 
local behavior of the approximating curves 7 e is chosen in a particular way, then both 
the local statistics and the global height fluctuations of S e can be determined. 
Here is the result on the local statistics. 

Theorem 1.2 Suppose that the gradient of h is not maximal at any point in U, that 
is, the normal vector to the surface has nonzero coordinates at every point. Under 
appropriate hypotheses on the local structure of the approximating curves %, the local 
statistics of T, e near a given point are given by the unique Gibbs measure of slope equal 
to the slope of h at that point. 
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Figure 2: The triangle and a scaled copy with vertices 0, 1, $. 

For the precise statement see Theorem 16.11 In particular the hypotheses on 7 e are 
explained in sections 12,61 and 13.21 

Suppose that (p a ,Pb,Pc) is a normal vector to the surface at a point. Recall that 
p a ,Pb,Pc > 0. If we rescale so that p a + Pb + Pc = F then one consequence of Theorem 
16. II is that the quantities p a ,Pb,Pc are the densities of the three orientations of lozenges 
near the corresponding point on the surface. 

1.2.3 Fluctuations 

The fluctuations are the image of the Gaussian free field under a certain diffeomorphism 
from the unit disk ID to U. To describe the fluctuations, we first describe the relevant 
conformal structure on U. It is a function of the normal to the graph of h, and is 
defined as follows. Let (p a ,Pb,Pc) be the normal to h, scaled so that p a + Pb + Pc = F 
Let 6 a = 7rp a ,9b = 7tpb,0 c = irp c . Let a,b,c be the edges of a Euclidean triangle 
with angles 9 a ,9b,0 c . Let z = —e~ l8c and w = —e lBb so that a + bz + cw = 0. See 
Figure [2J here the triangle on the left has edges a, bz, cw when these edges are oriented 
counterclockwise. We define $ = —cw/a. All of these quantities are functions on U, 
although a, b, c are only defined up to scale. Let x, y, z be the unit vectors in Pm in 
the directions of the projections of the standard basis vectors in M 3 . 

Theorem 1.3 (|llj) The function $ satisfies the complex Burgers equation 

$ £ + *$ fi = 0, (1) 

where &x,&y are directional derivatives of & in directions x,y respectively. 

The function : U — > C can be used to define a conformal structure on U, as 
follows. A function g : U — > C is defined to be analytic in this conformal structure 
on U if it satisfies g% + &gy = 0, where gx,gy are the directional derivatives of g in 
the directions x, y respectively. By the Alhfors-Bers theorem there is a diffeomorphism 
/ : U — ► B satisfying f x + &f y = 0; the conformal structure on U is the pull-back 
of the standard conformal structure on D under /. The conformal structure on U 
can be described by a Beltrami coefficient £ (see below) which in the current case is 
£ = ($ - e i7r / 3 )/($ - e~ i7r / 3 ). 
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In the special cases that 3> is constant (which correspond to the cases where 7 is 
contained in a plane), this means that the conformal structure is just a linear image of 
the standard conformal structure. For example, note that in the standard conformal 
structure on Pm, a function g is analytic if g% + e ln ^ 3 gy = 0. So the case <3? = e ln ^ 3 , 
which corresponds to the case a = b = c, gives the standard conformal structure (recall 
that the vectors x and y are 120° apart). 

Theorem 1.4 Suppose that the gradient ofh is not maximal at any point in U . Under 
the same hypotheses on 7 e as in Theorem the fluctuations of the unnormalized 
height function, -{h e — h), have a weak limit as e — > which is the Gaussian free field 
in the complex structure defined by <£, that is, the pull-back under the map f : U — > D 
above, of the Gaussian free field on the unit disk B. 

For the definition of the Gaussian free field see below. 

As mentioned above, we require that the normal to the graph of h be strictly 
inside the positive orthant, so that we have a positive lower bound on the values of 
PaiVbiPc, whereas the results of [3] and [UJ do not require this restriction. Indeed, in 
many of the simplest cases the surface So will have facets, which are regions on which 
(p a ,Pb,Pc) = (1) 0, 0), (0, 1,0) or (0, 0, 1). Our results do not apply to these situations. 

This work builds on work of [9[ \12\ ITT] . Previously Theorem 11.41 was proved for a 
the dimer model on Z, 2 in a special case which in our context corresponds to the wire 
frame 7 lying in the plane Pm, see [9]. In that case the conformal structure is the 
standard conformal structure on U. 

In the present case we still require special boundary conditions, which generalize 
the "Temperleyan" boundary conditions of [HE]. It remains an open question whether 
the result holds for all boundary conditions. The fluctuations in the presence of facets 
are also unknown, and the current techniques to not seem to immediately extend to 
this more general setting. 



1.3 The Gaussian free field 

The Gaussian free field X on B [T7] is a random object in the space of distributions on 
B, defined on smooth test functions as follows. For any smooth test function ip on B, 
J il){x)X(x)\dx\ 2 is a real Gaussian random variable of mean zero and variance given 

by 

'4){xi)'4){x2)G{xi,X2)\dxi\ 2 \dx2\ 2 



where the kernel G is the Dirichlet Green's function on I 

x\ - x 2 



G(x 1 ,x 2 ) = - — log 

Z7T 



1 - 



A similar definition holds (for the standard conformal structure) on any bounded 
domain in C, only the expression for the Green's function is different. 

An alternative description of the Gaussian free field is that it is the unique Gaus- 
sian process which satisfies E[X(xi)X(x2)] = G{x\,X2)- Higher moments of Gaussian 
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processes can always be written in terms of the moments of order 2; for the Gaussian 
free field we have 

E(X(xi)...X{x n )) = if n is odd, 

and 

E(X(x 1 )...X(x 2k )) = Yl G ( X v(l)' X v(2))--- G ( X *(2k-l),Xa(2k)) (2) 

pairings 

where the sum is over all (2k — 1)!! pairings of the indices. Any process whose moments 
satisfy ([2]) is the Gaussian free field [9]. 

1.4 Beltrami coefficient 

A conformal structure on U can be defined as an equivalence class of diffeomorphisms 
(f> : U — > D, where mappings (f>i,4>2 are equivalent if the composition <p\ o cj)^ 1 is a 
conformal self-map of D. The Beltrami differential of <f> is defined by the 

formula 

£(z)— - ti— 

dz 4> z dz 

The Beltrami differential is invariant under post-composition of <j> with a conformal 
map, so it is a function only of the conformal structure (and in fact defines the conformal 
structure as well). It is not hard to show that \£,(z)\ < 1; note that = if and 
only if the map is conformal. The Ahlfors-Bers uniformization theorem [I] says that 
any smooth function (even any measurable function) £(z) satisfying |£(z)| < 1 defines 
a conformal structure. 

1.5 Examples 

The simplest case is when the wire frame 7 is contained in a plane {(x, y, z) G M 3 | p a x+ 
PbV + Pc z = const}. In this case the limit surface £0 is linear. The normal (p a ,Pb,Pc) 
is constant, and the conformal structure is a linear image of the standard conformal 
structure. That is, the map / : U — > B is a linear map L composed with a conformal 
map. 

For a more interesting case, consider the boxed plane partition (BPP) shown 
in Figure [3l which is a random lozenge tiling of a regular hexagon. In [1] it was 
shown that, for a random tiling of the hexagon, the asymptotic height function h is 
linear outside of the inscribed circle and analytic inside (with an explicit but somewhat 
complicated formula). Although our theorem does not apply to this case because of 
the facets outside the inscribed circle, if we choose boundary conditions inside the 
inscribed circle, and boundary values equal to the graph of the function there, our 
results apply. Suppose that the hexagon has sides of length 1, so that the inscribed 
circle has radius y/3/2. Let U be a disk of radius r < y/3/2 concentric with it. Suppose 
that the normalized height function on the boundary of U is chosen to agree with 
the asymptotic height function of the corresponding region in the BPP, so that the 
asymptotic height function of U equals the asymptotic height function of the BPP 
restricted to U. Then the fluctuations on U can be computed using Theorem 17.11 
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Figure 3: Boxed plane partition. 
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In fact in this setting, the conformal structure can be explicitly computed: take the 
standard conformal structure on a hemisphere in M 3 , and project it orthogonally onto 
the plane containing its equator. Identifying the equator with the inscribed circle in 
the BPP gives the relevant conformal structure in U. Remarkably, in this example the 
Beltrami coefficient is rotationally invariant, even though h itself is not. See section [HJ 
We conjecture that the fluctuations for the BPP are given by the limit r — > v3/2 
of this construction (it is known [3] that fluctuations in the "frozen" regions outside 
the circle are exponentially small in 1/e). 



1.6 Proof outline 

The fundamental tool in the study of the dimer model is the Kasteleyn matrix (de- 
fined below). Minors of the inverse Kasteleyn matrix compute edge correlations in 
the model. The main goal of the paper is to obtain an asymptotic expansion of the 
inverse Kasteleyn matrix. This is complicated by the fact that it grows exponentially 
in the distance between vertices (except in the special case when the boundary height 
function is horizontal). However by pre- and post-composition with an appropriate di- 
agonal matrix, we can remove the exponential growth and relate K^ 1 to the standard 
Green's function with Dirichlet boundary conditions on a related graph St- 
Here is a sketch of the main ideas. 

1. We construct a discrete version of the map / of Theorem ll.41 For each e we define 
a directed graph St embedded in the upper half plane H, and a geometric map 
(f) from U e to St, such that the Laplacian on St is related (via the construction 
of |13|. [14] ) with the Kasteleyn matrix on U e . The existence of such a graph St 
follows from p3j. This is done in section [3.11 in the "constant slope" case and 
section 14.31 in the general case. 

2. Standard techniques for discrete harmonic functions yield an asymptotic expan- 
sion for the Green's function on St- This is done in section [2.6.41 

3. The asymptotic expansion of the inverse Kasteleyn matrix on U e is obtained from 
the derivative of the Green's function on Stj pulled back under the mapping <p. 
See sections 12.6.41 and [6j 

4. Asymptotic expansions of the moments of the height fluctuations are computed 
via integrals of the asymptotic inverse Kasteleyn matrix. These moments are the 
moments of the Gaussian free field on EI pulled back under (f>. 

Acknowledgments. Many ideas in this paper were inspired by conversations with 
Henry Cohn, Jim Propp, Jean-Rene Geoffroy, Scott Sheffield, Beatrice deTiliere, Cedric 
Boutillier, and Andrei Okounkov. We thanks the referees for useful comments. This 
paper was partially completed while the author was visiting Princeton University. 
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2 Definitions 



2.1 Graphs 

Let 7riii be the orthogonal projection of M 3 onto Pm. Let x, y, z be the 7rm-projections 
of the unit basis vectors. Define e% = |(x — z), &2 = — x), e% = — y), so that 
x = e\ — e2, y = &2 — e3,z = e% — e%. Let be the honeycomb lattice in Pm: vertices 
of K are LU(I + ei), where L is the lattice L = Z(ei — e2) + Z(e2 — 63) = Zx + Zy, 
and edges connect nearest neighbors. Vertices in L are colored white, those in L + e\ 
are black. See Figured! 



Z X 



Figure 4: The honeycomb graph. 



2.1.1 Dual graph 

Let U be a Jordan domain in Pm with smooth boundary. In e!K, take a simple closed 
polygonal path with approximates dU in a reasonable way, for example the polygonal 
curve is locally monotone in the same direction as the curve dU. Let 9 be the subgraph 
of e'K bounded by this polygonal path. We define a special kind of dual graph 9* as 
follows. Let e'K* be the usual planar dual of e'K. For each white vertex of 9 take the 
corresponding triangular face of eK*; the union of the edges forming these triangles, 
along with the corresponding vertices, forms 9*- In other words, 9* has a face for each 
white vertex of 9, as well as for black vertices which have all three neighbors in 9- See 
Figure El 

Throughout the paper the graph 9 and its related graphs will be scaled by a factor 
e over the corresponding graphs K, and so the 9 graphs have edge lengths of order e, 
and the graph K and its related graphs have edge lengths of order 1. 

2.1.2 Forms 

For an edge in 9 joining vertices b and w, we denote by (bw)* the dual edge in 9*, 
which we orient at +90° from the edge bw (when this edge is oriented from b to w). 

A 1-form u on a graph is a function on directed edges which is antisymmetric with 
respect to reversing the orientation: u(viV2) = —uj(v2Vi). A 1-form is also called a 
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Figure 5: The "dual" graph S* (solid lines) of the graph 9 of Figure HI 



flow. If the graph is planar one can similarly define a 1-form on the dual graph. If u 
is a 1-form, w* is the dual 1-form, defined by w* {(v\V2)*) '■= w(v\V2)- 

On a planar graph with a 1-form w, dw is a function on oriented faces defined by 
dto(f) = X^ e u; ( e ) where the sum is over the edges on a path around the face going 
counterclockwise. This is also known as the curl of the flow w. The form dw* is a 
function on vertices (faces of the dual graph), defined by dw*(v) = Yl v '~ v u; (. v > v ')- ^ n 
other words it is the divergence of the flow u. 

A 1-form w is closed if dw = 0, that is, the sum of ui along any cycle is zero (in the 
language of flows, the flow has zero curl). If is exact if oj = df for some function / on 
the vertices, that is u>(viV2) = f(vi) — f{ v 2)- A 1-form is co-closed if its dual form is 
closed. The corresponding flow is divergence- free. 

If dw* = 0, the integral of w* between two faces of 9 (i-e. on a path in the dual 
graph) is the flux, or total flow, between those faces. 

2.2 Heights and asymptotics 

The unnormalized height function, or just height function, of a tiling is the integer- 
valued function on the vertices of the lozenges (faces of 9) which is the sum of the 
coordinates of the corresponding point in Z 3 . It changes by ±1 along each edge of a 
tile. When we scale the lattice by e, so that we are discussing surfaces in eZ 3 , the height 
function is defined as 1/e times the coordinate sum, so that it is still integer- valued. 
The normalized height function is e times the height function, and is the function 
which, when scaled by y/3, has graph which is the surface in eZ 3 . 

Let u : dU ->lbea continuous function with the property that u can be extended 
to a Lipschitz function u on the interior of U having the property that the normal to 
the graph of u has nonnegative coordinates, that is, the normal points into the positive 
orthant M 3 . In other words, the graph of u is a wire frame 7 of the type discussed 
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before. 

Let h : U —> R be the asymptotic height function with boundary values u, from 
Theorem 11.11 It is smooth assuming the hypothesis of Theorem 11.41 

Let v = (j> a ,Pb,Pc) be the normal vector to the graph of h, scaled so that p a +Pb + 
p c = 1. The directional derivatives of h in the directions x,y,z are respectively 

3p c - 1, 3p a - 1, 3p b - 1. (3) 

2.3 Measures and gauge equivalence 

We let fjL = /i(S) be the uniform measure on dimer configurations on a finite graph 9- 

If edges of 9 are given positive real weights, we can define a new probability measure, 
the Boltzmann measure, giving a configuration a probability proportional to the 
product of its edge weights. 

Certain edge-weight functions lead to the same Boltzmann measure: in particular 
if we multiply by a constant the weights of all the edges in 9 having a fixed vertex, 
the Boltzmann measure does not change, since exactly one of these weights is used in 
every configuration. More generally, two weight functions u\ , v<i are said to be gauge 
equivalent if v\jv2 is a product of such operations, that is, if there are functions F\ 
on white vertices and F2 on black vertices so that for each edge wb, i/i(wb)/z/2(wb) = 
i ? i(w)i ? 2(b). Gauge equivalent weights define the same Boltzmann measure. 

It is not hard to show that for planar graphs, two edge-weight functions are gauge 
equivalent if and only if they have the same face weights, where the weight of a face 
is defined to be the alternating product of the edge weights around the face (that is, 
the first, divided by the second, times the third, and so on), see e.g. |12j . 

In this paper we will only consider weights which are gauge equivalent to constant 
weights (or nearly so), so the Boltzmann measure will always be (nearly) the uniform 
measure. 

2.4 Kasteleyn matrices 

Kasteleyn showed that one can count dimer configurations on planar graph with the 
determinant of the certain matrix, the "Kasteleyn matrix" |6j. In the current case, 
when the underlying graph is part of the honeycomb graph, the Kasteleyn matrix K 
is just the adjacency matrix from white vertices to black vertices. 

For more general bipartite planar graphs, and when the edges have weights, the 
matrix is a signed, weighted version of the adjacency matrix [16], whose determinant 
is the sum of the weights of dimer coverings. Each entry K(w, b) is a complex number 
with modulus given by the corresponding edge weight (or zero if the vertices are not 
adjacent), and an argument which must be chosen in such a way that around each face 
the alternating product of the entries (the first, divided by the second, times the third, 
and so on) is positive if the face has 2 mod 4 edges and negative if the face has mod 4 
edges (since we are assuming the graph is bipartite, each face has an even number of 
edges. For nonbipartite graphs, a more complicated condition is necessary). 
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The Kasteleyn matrix is unique up to gauge transformations, which consist of pre- 
and post-multiplication by diagonal matrices (with, in general, complex entries). If the 
weights are real then we can choose a gauge in which K is real, although in certain 
cases it is convenient to allow complex numbers (we will below). 

Probabilities of individual edges occurring in a random tiling can likewise be com- 
puted using the minors of the inverse Kasteleyn matrix: 

Theorem 2.1 ([7]) The probability of edges {(t>i, wi), . . . , (bfe, w&)} occurring in a 
random dimer covering is 



On an infinite graph K is defined similarly but K is not unique in general. This 
is related to the fact that there are potentially many different measures which could 
be obtained as limits of Boltzmann measures on sequences of finite graphs filling out 
the infinite graph. The edge probabilities for these measures can all be described as in 
the theorem above, but where the matrix "X -1 " n0 w depends on the measure; see the 
next section for examples. 

2.5 Measures in infinite volume 

On the infinite honeycomb graph "K there is a two-parameter family of natural translation- 
invariant and ergodic probability measures on dimer configurations, which restrict to 
the uniform measure on finite regions (i.e. when conditioned on the complement of 
the finite region: we say they are conditionally uniform). Such measures are also 
known as ergodic Gibbs measures. They are classified in the following theorem due 
to Sheffield. 

Theorem 2.2 ( [18J ) For each u = (p a iPb,Pc) with p a ,Pb,Pc > and scaled so that 
Pa + Pb + Pc = 1 there is a unique translation-invariant ergodic Gibbs measure fj, v on 
the set of dimer coverings of "K, for which the height function has average normal v. 
This measure can be obtained as the limit as n — > oo of the uniform measure on the 
set of those dimer coverings of "K n = "K/nL whose proportion of dimers in the three 
orientations is (p a : pb ■ p c )> up to errors tending to zero as n —* oo. Moreover every 
ergodic Gibbs measure on "K is of the above type for some v. 

The unicity in the above statement is a deep and important result. 

Associated to \i v is an infinite matrix, the inverse Kasteleyn matrix of fj, u , 
K^ 1 = (K~ l (b,w)) whose rows index the black vertices and columns index the white 
vertices, and whose minors give local statistics for /i u , just as in Theorem 12.11 From 
|12j there is an explicit formula for K„ : let w = m\x + n\y and b = e\ + m2X + niy = 
w + e\ + rax + ny where m = — mi,n = n-i — n\. Then 





(4) 
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where a, b, c, z, w are as defined in section f)1.2.3j) and 



K-\b w) = I H m+n w~ x n dz 1 dw 1 

abc ' (2iri) 2 J\ z \=\ w \ = i a + bz\ + cw\ z\ w\ 

_ I to (£^) +0 (_i_). (6) 

7T \ cwm + on / \m z + n z J 

This formula for if"T and its asymptotics were derived in |12| : they are obtained 
from the limit n — ► oo of the inverse Kasteleyn matrix on the torus "K/nL with edge 
weights a, b, c according to direction. It is not hard to check from ([5]) that KK^ 1 = Id 

From ([3]), the matrix K^ 1 is just a gauge transformation of K~ bc , that is, obtained 
by pre- and post-composing with diagonal matrices. 

Defining F(w) = (bz/a) mi (cw/bz) ni and F(b) = a(bz/a)~ m2 (cw/bz)~ n2 we can 
write, using Q, 



\ cwm + an J m z + n z 



W)|P( - ; , ). (7) 



We'll use this function T below. 

As a sample calculation, the /i^-probability of a single horizontal edge, from w = 
to b = ei , being present in a random dimer covering is (see Theorem 12. ip 



^(b,w)=a^(b,w) = ^jf 



1 dz\ dw\ 9 a 

j2 a + 6zi + cw\ Z\ w\ it 



where 9 a is, as before, the angle opposite side a in a triangle with sides a, b, c. This is 
consistent with ([3|). 



2.6 T- graphs 
2.6.1 Definition 

T-graphs were defined and studied in [H]- A pairwise disjoint collection Li, L2, . . . ,L n 
of open line segments in M 2 forms a T-graph in M 2 if U™ =1 Lj is connected and 
contains all of its limit points except for some finite set R = {n, . . . , r m }, where each 
rj lies on the boundary of the infinite component of M? minus the closure of U™ =1 Lj. 
See Figure [6] for an example where the outer boundary is a polygon. Elements in R 
are called root vertices and are labeled in cyclic order; the Lj are called complete 
edges. We only consider the case that the outer boundary of the T-graph is a simple 
polygon, and the root vertices are the convex corners of this polygon. (An example 
where the outer boundary is not a polygon is a "T" formed from two edges, one ending 
in the interior of the other.) 

Associated to a T-graph is a Markov chain St, whose vertices are the points which 
are endpoints of some Lj. Each non-root vertex is in the interior of a unique Lj (because 
the Lj are disjoint); there is a transition from that vertex to its adjacent vertices along 
Lj, and the transition probabilities are proportional to the inverses of the Euclidean 
distances. Root vertices are sinks of the Markov chain. See Figure 
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Figure 7: The Markov chain associated to the T-graph of Figure El Note that root vertices 
are sinks of the Markov chain. 
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Note that the coordinate functions on St are harmonic functions on St \ R- More 
generally, any function / on Sr which is harmonic on St\R (we refer to such functions 
as harmonic functions on St) has the property that it is linear along edges, that is, if 
v\ , V2 , v 3 are vertices on the same complete edge then 

f(vi) - f_M = /Mi(M , (8) 

Vl - V 2 V2 — V3 

2.6.2 Associated dimer graph and Kasteleyn matrix 

Associated to a T-graph is a weighted bipartite planar graph Sd constructed as follows, 
see Figure EJ Black vertices of Sd are the Lj. White vertices are the bounded comple- 
mentary regions, as well as one white vertex for each boundary path joining consecutive 
root vertices rj and r^+i (but not for the path from r m to r\). The complementary re- 
gions are called faces; the paths between adjacent root vertices are called outer faces. 
Edges connect the Li to each face it borders along a positive- length subsegment. The 
edge weights are equal to the Euclidean length of the bounding segment. 

To Sd there is a canonically associated Kasteleyn matrix of Sd : this is the n x n 
matrix Kg D = (Kg D (w, b)) with rows indexing the white vertices and columns indexing 
the black vertices of Sd- We have ifg D (w, b) = if w and b are not adjacent, and 
otherwise Kg D (w, b) is the complex number equal to the edge vector corresponding to 
the edge of the region w along complete edge b (taken in the counterclockwise direction 
around w). In particular |JCg (w, b)| is the length of the corresponding edge of w. 

Lemma 2.3 Kg D is a Kasteleyn matrix for S d , that is, the alternating product of the 
matrix entries for edges around a bounded face is positive real or negative real according 
to whether the face has 2 mod 4 or mod 4 edges, respectively. 

By alternating product we mean the first, divided by the second, times the third, 

etc. 

Proof: Let / be a bounded face of Sd (we mean not one of the outer faces). It 
corresponds to a meeting point of two or more complete edges; this meeting point is 
in the interior of exactly one of these complete edges, L. See Figure El In 9d> for 
each other black vertex on that face the two edges of the T-graph to neighboring white 
vertices have opposite orientations. The two edges parallel to L (horizontal in the 
figure) have the same orientation, so their ratio is positive. This implies the result. □ 

Although we won't need this fact, in [14] it is shown that the set in-directed spanning 
forests of Sr (rooted at the root vertices and weighted by the product of the transition 
probabilities) is in measure-preserving (up to a global constant) bijection with the set 
of dimer coverings of Sd ■ 

2.6.3 Harmonic functions and discrete analytic functions 

To a harmonic function / on a T-graph St we associate a derivative df which is a 
function on black vertices of Sd as follows. Let v\ and V2 be two distinct points on 
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Figure 8: A face of Sd (dotted). 



complete edge b, considered as complex numbers. We define 

dm = /(^WM , (9) 

Since / is linear along any complete edge (equation ([SD), df is independent of the choice 
of v\ and V2- 

Lemma 2.4 /// is harmonic on a T -graph $t and K = Kg D is the associated Kaste- 
leny matrix, then ^2 beB K(w, b)df(b) = for any interior white vertex w. 

Proof: Let bx, . . . , b& be the neighbors of w in cyclic order. To each neighbor bj is 
associated a segment of a complete edge Lj. Let Vi and Vi + \ be the endpoints of that 
segment, and Wi and w[ be the endpoints of Lj. Then 

f( Vl+1 ) - fJVj) _ /K) - f{Wj) 
V i+ i - Vi W • - Wi 

since the harmonic function is linear along Lj. In particular 

k Sd ( w , b t )df(b t ) = (v t+1 - = f(Vt+l) _ 

tfj Wi 

Summing over i (with cyclic indices) yields the result. □ 

Note that at a boundary white vertex w, Kc ID (w, b)d/(b) is the difference in 
/-values at the adjacent root vertices of St- 

We will refer to a function g on black vertices of Sd satisfying X^be-B K-Sd ( w ' b)g(b) = 
for all interior white vertices w as a discrete analytic function. 

The construction in the above lemma can be reversed, starting from a discrete an- 
alytic function df (on black vertices of Sd) an d integrating to get a harmonic function 
/ on St : define / arbitrarily at a vertex of St and then extend to neighboring ver- 
tices (on a same complete edge) using Q. The extension is well-defined by discrete 
analyticity. 
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Figure 9: Conjugate Green's function example. Here the sides of the triangle are bisected in 
ratio 2 : 3 and interior complete edges 1:1. 

2.6.4 Green's function and 

We can relate to the conjugate Green's function on St using the construction of 
the previous section, as follows. 

Let w be an interior face of St> and I a path from a point in w to the outer boundary 
of St which misses all the vertices of St- For vertices v of St, define the conjugate 
Green's function G*(w,v) to be the expected algebraic number of crossings of £ 
by the random walk started at v and stopped at the boundary. This is the unique 
function with zero boundary values which is harmonic everywhere except for a jump 
discontinuity of —1 across I when going counterclockwise around w. (If there were 
two such functions, their difference would be harmonic everywhere with zero boundary 
values. ) See Figure [9] for an example. 

Let Kg^(b, w) be the discrete analytic function of b defined from G*(w, v) as in the 
previous section; on an edge which crosses £, define K^(b,w) using two points on b 
on the same side of I. This function clearly satisfies K^ d K^ d = I by Lemma 12.41 and 
therefore is independent of the choice of I. 

Note that also has the following probabilistic interpretation: take two particles, 
started simultaneously at two different points v±,V2 of the same complete edge, and 
couple their random walks so that they start independently, take simultaneous steps, 
and when they meet they stick together for all future times. Then the difference in 
their winding numbers around w is determined by their crossings of I before they meet. 
That is, Kg (b,w)(«i — V2) is the expected difference in crossings before the particles 
meet or until they hit the boundary, whichever comes first. 
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3 Constant-slope case 



In this section we compute the asymptotic expansion of Kg* (Theorem 13. 7p in the 
special case is when 7 is planar. In this case the normalized asymptotic height function 
h is linear, and its normal v is constant. This case already contains most of the 
complexity of the general case, which is treated in section 01 

Let v = (p a ,Pb,Pc) be the normal to h, scaled as usual so that p a + Pb + Pc = 1- 
The angles 6 a , 6b, 9 C are constant, and we choose a, b, c as before to be constant as well. 
Define a function 

F(w) = {bz/a) m {cw/bz) n 
at a white vertex w = mi + ny and 

F(b) =a(bz/a)- m (cw/bz)- n 

at a black vertex b = ei + mx+ny. These functions are defined on all of the honeycomb 
graph "K. Let Kj^ be the adjacency matrix of "K (which, as we mentioned earlier, is a 
Kasteleyn matrix for 'K). 

Lemma 3.1 We have 

^^(w,b)F(b) =0 = ^F(w)K M (w,b) (10) 

b w 

where Kj£ is the adjacency matrix of Ji. 

Proof: This follows from the equation a + bz + cw = 0. □ 



3.1 T-graph construction 

Define a 1-form $7 on edges of !K by 

fi(wb) = -O(bw) = 2Re(F(w))F(b). (11) 

By (JTOj) the dual form $7* (defined by fT((wb)*) = f2(wb)) is closed (the integral 
around any closed cycle is zero) and therefore fT = for a complex-valued function ^ 
on "K* . Here "K* , the dual of the honeycomb, is the graph of the equilateral triangulation 
of the plane. Extend ^ linearly over the edges of "K* . This defines a mapping from !K* 
to C with the property that the images of the white faces are triangles similar to the 
a, b, c-triangle (via orientation-preserving similarities), and the images of black faces 
are segments. This follows immediately from the definitions: if bi, b2, b3 are the three 
neighbors of a white vertex w of JC, the edges wbj have values 2Re(i ? (w))i ? (bj) which 
are proportional to F(b\) : F(b2) ■ F(bs), which in turn are proportional to a : bz : cw 
by (jlOp . If wi, W2, W3 are the three neighbors of a black vertex then the corresponding 
edge values are 2Re(F(wj))i ? (b) which are proportional to Re(F(wi)) : Re(F(\«2)) : 
Re(-F(w3)), that is, they all have the same slope and sum to zero. 
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Figure 10: The \l/-image of %* is a T-graph covering IR 2 



It is not hard to see that the images of the white triangles are in fact non-overlapping 
(see |14] , Section 5 for the proof, or look at Figure [TO]) . It may be that Rei ? (w) = for 
some w; in this case choose a generic modulus-1 complex number A and replace F(w) 
by A-F(w) and F(b) by AF(b). So we can assume that each white triangle is similar to 
the a, b, c-triangle. In fact, this operation will be important later; note that by varying 
A the size of an individual triangle varies; by an appropriate choice we can make any 
particular triangle have maximal size (side lengths a,b,c). 

Lemma 3.2 The mapping ^ is almost linear, that is, it is a linear map 4>(m, n) = 
cwm + an plus a bounded function. 

Proof: Consider for example a vertical column of horizontal edges {wibi, . . . , w^b/^} 
of "K connecting a face f\ to face fk+i- We have 
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^0*((w i b i )*) = ^fi^bi) 

i=i i=i 

k 

= ^( J P(w i ) + FK))F(b i ) 

i=l 

k 

= fca + ^F(^)F(bi) (12) 



i=l 
k 



ka + a 



z 2l w- 21 



i=l 

= 4>(fk+l) ~ 4>{fl) +OSC, 

where osc is oscillating and in fact O(l) independently of k (by hypothesis z/w M). 
In the other two lattice directions the linear part of ^ is again 0, so that ^ is almost 
linear in all directions. □ 

This lemma shows that the image of !H* under is an infinite T-graph "Kt covering 
all of M 2 . The images of the black triangles are the complete edges and have lengths 
O(l). 

If we insert a A as above and let A vary over the unit circle, one sees all possible 
local structures of the T-graph, that is, the geometry of the T-graph "Kt = !Kt(A) in 
a neighborhood of a triangle ^(w) only depends up to homothety on the argument of 
AF(w). 

Recall that S* is a subgraph of e'K* approximating U. We can restrict ^ to -S* 
thought of as a subgraph of and then multiply its image by e. Thus we get a finite 
sub-T-graph S T of e%r- Let ^ e = e o $ o ± so that ^ e acts on S*. 

The union of the ^-images of the white triangles in S* forms a polygon P. Define 
the "dimer" graphs !Hd associated to "Kt, and Sd associated to St as in section 12.6.21 
See Figure [6] for the T-graph arising from the graph S* of Figure Note that S d 
contains S (defined in section [2. l.lj) but has extra white vertices along the boundaries. 
These extra vertices make the height function on 3d approximate the desired linear 
function (whose graph has normal u), see the next section. 

From (1111) we have 



Lemma 3.3 Edge weights of^Hn and Sr> are gauge equivalent to constant edge weights. 
Indeed, for w not on the boundary of Sd we have 

K Sd (w, b) = 2eRe(T(w))F(b)^ s (w, b) (13) 

and similarly for all w, 

Kk d = 2Re(T(w))T(b)^(w, b). (14) 
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Asymptotics of K q 1 are described below in section 13.41 For K^c , from ([7]) and 
141) we have 



= 2,Re(F( 1 w))F(b) Im - «*) ) + 0> — °*> 



1 

|0(b) ~0(w)P 



Here -K"-^ is the inverse constructed from the conjugate Green's function on the T- 
graph Ji T . We use the fact that 2Re(F(w))F(b)iC^ (b, w) coincides with K" 1 of © 
since both satisfy the equation that dK -1 equals the conjugate Green's function. 



3.2 Boundary behavior 

Recall that from our region U we constructed a graph 9 (section |2.1.ip . From the 
normalized height function u on dU (which is the restriction of a linear function to 
dU) we constructed the T-graph St and then the dimer graph Sd- 

In this section we show that the normalized boundary height function of Sr> when 
3d is chosen as above approximates h. This is proved in a roundabout way: we first 
show that the boundary height does not depend (up to local fluctuations) on the exact 
choice of boundary conditions, as long as we construct St from S* as in section 12.11 
Then we compute the height change for "simple" boundary conditions. 



3.2.1 Flows and dimer configurations 

Any dimer configuration m defines a flow (or 1-form) [m] with divergence 1 at each 
white vertex and divergence —1 at each black vertex: just flow by 1 along each edge in 
m and on the other edges. 

The set £l\ C [0, 1] E of unit white-to-black flows (i.e. flows with divergence 1 at 
each white vertex and —1 at each black vertex) and with capacity 1 on each edge is 
a convex polytope whose vertices are the dimer configurations [15] . On the graph "K 
define the flow 0^/3 6 f2i to be the flow with value 1/3 on each edge wb from w to b. 
Up to a factor 1/3, this flow can be used to define the height function, in the sense that 
for any dimer configuration m, [m] — is a divergence- free flow and the integral of 
its dual (which is closed) is 1/3 times the height function of m. That is, the height 
difference between two points is three times the flux of [m] — W1/3 between those points. 
This is easy to see: across an edge which contains a dimer the height changes by ±2 
(depending on the orientation) ; if the edge does not contain a dimer the height changes 
by =Fl, for so that the height change can be written i(3% — 1) where x 1S the indicator 
function of a dimer on that edge. 

For a finite subgraph of "K the boundary height function can be obtained by inte- 
grating around the boundary (3 times) the dual of [m] — W1/3 for some m. 



3.2.2 Canonical flow 

On !Hd there is a canonical flow cu £ f^i defined as illustrated in Figure [TO let V1V2 be 
the two vertices of St along b adjacent to w. The flow from w to b is l/(27r) times the 
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sum of the two angles that the complete edges through v\ and V2 make with b. One 
or both of these angles may be zero. 

Note that the total flow out of w is 1. For an edge b, the total flow into b is also 1 
as illustrated on the right in Figure PTT1 




Figure 11: Defining the canonical flow (divide the angles by 2tt). 

Now Sd is a subgraph of CKd with extra white vertices around its boundary. The 
canonical flow on IKd restricts to a flow on Sd except on edges connecting to these 
extra white vertices (we define the canonical flow there to be zero). This flow has 
divergence 1,-1 at white/black vertices, except at the black vertices of Sd connected 
to the boundary white vertices, and the boundary white vertices themselves. If m is a 
dimer covering of 3d > the flow [to] — uj is now a divergence- free flow on Sd except at 
these black and white boundary vertices. 

Lemma 3.4 Along the boundary of Sd the divergence of [m] — uj for any dimer con- 
figuration to is the turning angle of the boundary of P. 

Proof: Consider a complete edge L corresponding to black vertex b. The canonical 
flow into b has a contribution from the two endpoints of L. The flow [m] can be 
considered to contribute —1/2 for each endpoint. 

Recall that each endpoint of L ends in the interior of another complete edge or at 
a convex vertex of the polygon P. If an endpoint of L ends in the interior of another 
complete edge, and there are white faces adjacent to the two sides of this endpoint (that 
is, the endpoint is not a concave vertex of P) then the contribution of the canonical 
flow is also —1/2, so the contribution of [to] — uj is zero. 

Suppose the endpoint is at a concave vertex v of P with exterior angle 9 < ir. The 
contribution from v to the flow of [m] — uj into L is— 1/2 -\- J- . The other complete 
edge at v does not end at v and so has no contribution. This quantity is 1/27T times 
the turning angle of the boundary. 

Suppose the endpoint is at a convex vertex v' of P of interior angle 6 < ir. There 
is some complete edge V of "Kt, containing v' in its interior, which is not in St- The 
sum of the contributions of [m] — uj for the endpoints of the two complete edges L%, L2 
of St meeting at v' is —1/2 — 9/2-k. This is the 1/27T times the turning angle at the 
convex vertex, minus 1. 
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The contribution from [m] from the boundary white vertices is 1 per white boundary 
vertex, that is, one per convex vertex of P. □ 

This lemma proves in particular that the divergence of [m] — lo is bounded for any 
[m]. 

3.2.3 Boundary height 

Recall that the height function of a dimer covering m can be defined as the flux of 
u>]/3 — [m]. In particular, the flux of uu 3 — u defines the boundary height function 
up to O(l) (the turning of the boundary), since the flux of [m] — uj is the boundary 
turning angle. 

The flux of LJ1/3 — uj between two faces can be computed along any path, and in 
fact because both uj 1 / 3 and uj are locally defined from "Hd, we see that the flux does 
not depend on the choice of the nearby boundary. Let us compute this flux and show 
that it is linear along lattice directions, and therefore linear everywhere in "Kd- 

Take a vertical column of horizontal edges in 9 , and let us compute uj on this set of 
edges. The ^ image of the dual of this column is a polygonal curve r] whose jth edge 
is (using (|12() ) a constant times 1 + (^) 2j - The image of the triangles in the vertical 
column of S* is as shown in Figure [T2l 



Figure 12: A column of triangles from St- 

The jth edge is part of a complete edge corresponding to the jth black vertex in 
the column. By the argument of the previous section, the flux is equal to the number 
of convex corners of this polygonal curve rj, that is, corners where the curve turns left. 
The curve rj has a convex corner when 

V 1 + (z/w) 2 i J ~ 

that is, using z/w = —e tSa , when 2j9 a £ [vr,7r + 29 a ]. Assuming that 9 a is irrational, 
this happens with frequency so the flux of (Ju$ - co along a column of length n 
is — |), and the average flux per edge is p a — 1/3. Therefore the average height 
change per horizontal edge is 3p a — 1. If 9 a is rational, a continuity argument shows 
that 3p a — 1 is still the average height change per horizontal edge. A similar result holds 
in the other directions, and (|3|) show that the average height has normal (p a ,Pb,Pc) as 
desired. 
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3.3 Continuous and discrete harmonic functions 



To understand the asymptotic expansion of G* , the conjugate Green's function on 9t> 
from which we can get K -1 , we need two ingredients. We need to understand the 
conjugate Green's function on JCt, and also the relation between continuous harmonic 
functions on domains in M? and discrete harmonic functions on (domains in) e'Kx- 



3.3.1 Discrete and continuous Green's functions 

On "Kti the discrete conjugate Green's function G*(w,v), for w G IK and v £ JCy, can 
be obtained from integrating the exact formula for K^ 1 given in as discussed in 
Section l2iT31 

As we shall see, this formula differs even in its leading term from the continuous 
conjugate Green's function, due to the singularity at the diagonal. Basically, because 
our Markov chain is directed, a long random walk can have a nonzero expected winding 
number around the origin. This causes the conjugate Green's function, which measures 
this winding, to have a component of Re log v as well as a part Imlogv. 

The continuous conjugate Green's function on the whole plane is 

g*(vi,v 2 ) = -f-Imlog(u 2 - vi). 

(We use G* to denote the discrete conjugate Green's function and g* the continuous 
version.) 

To compute the discrete conjugate Green's function G* , if we simply restrict g* 
to the vertices of eJCr, it will be very nearly harmonic as a function of 1)2 for large 
\ v 2 ~ v i\, but the discrete Laplacian of g* at vertices v 2 near v\ (within O(e) of v\) will 
be of constant order in general. We can correct for the non-harmonicity at a vertex 
v' by adding an appropriate multiple of the actual (non-conjugate) discrete Green's 
function G(y,V2). The large scale behavior of this correction term is a constant times 
Relog(t>2 — v'). So we can expect the long-range behavior of the discrete conjugate 
Green's function G*(w, v), for |w — v\ large, to be equal to g*(w, v) plus a sum of terms 
involving the real part Relog(?; — v') for v's within 0(e) of <fr(w). These extra terms 
sum to a function of the form clog|v — + es(v) + 0(e 2 ), where s is a smooth 

function and c is a real constant, both c and s depending on the local structure of 'Kt 
near </>(w). 

This form of G* can be seen explicitly, of course, if we integrate the exact formula 
(fT5l). We have 



Lemma 3.5 The discrete conjugate Green's function G* on the plane e!K-r is asymp- 
totically 

G*(w,v) = i- (lmlog(« - 0(w)) + ^[^ Relogfo - 0(w))) + 

+ es( W ,v) + 0(e 2 ) + 0(- (16) 

\v — <p(w)\ 
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where s(w, v) is a smooth function of v. 



) (17) 



Proof: From the argument of the previous paragraph, it suffices to compute the 
constant in front of the Re \og(v — 4>(w)) term. This can be computed by differentiating 
the above formula with respect to v and comparing with formula (|15p for K~ l . The 
differential of es(w, v) is 0(e). Let v±,V2 be two vertices of complete edge b, coming 
from adjacent faces of "K, adjacent across an edge bwi (so that v i — v 2 = f2*(bwi) = 
2eRe(i ? (wi))i ? (b)). We have for w far from wi 



K^(b,w) 



G*(w, Vl )-G*(w,v 2 ) 



Vl - V 2 

V1-V2 \ 

+ 0(e) 



J_Tm ( F ( w ) "1~«2 
2vr 1 ^Rc(F(w)) 0(b)-0(w) 



2eRe(F(wi))F(b) 

1 / F(w)F(b) \ 

2vrRe(F(w))F(b) ™ \</>(b) - <£(w) J + [€) 



□ 



Note that in fact for any complex number A of modulus 1, we get a discrete conjugate 
Green's function G\ on the graph eJi/r(A) (from Section fe.lj) with similar asymptotics. 

3.3.2 Smooth functions 

Constant and linear functions on M 2 , when restricted to eJfy, are exactly harmonic. 
More generally, if we take a continuous harmonic function / on M. 2 and evaluate it on 
the vertices of e'Kx, the result will be close to a discrete harmonic function f 6 , in the 
sense that the discrete Laplacian will be 0(e 2 ): if v is a vertex of e'Kx and t>i,W2 are 
its (forward) neighbors located at v\ = v — ed\e lB and v 2 = v + ecfoe* 6 ' then the Taylor 
expansion of / about v yields 

Af(v) = f(v) - -j-^-rHv - ed ie w ) - + ed 2 e ie ) = 0{e 2 ). 

This situation is not as good as in the (more standard) case of a graph like el? , where 
if we evaluate a continuous harmonic function on the vertices, the Laplacian of the 
resulting discrete function is 0(e 4 ): 

f(v) ~ \ (f(v + 6) + f(v + ei) + f(v -e) + f(v - «)) = 0(e 4 ). 

In the present case the principal error is due to the second derivatives of /. To 
get an error smaller than 0(e 2 ), we need to add to / a term which cancels out the 
e 2 error. We can add to / a function which is e 2 times a bounded function / 2 whose 
value at a point v depends only on the local structure of e%r near v and on the second 
derivatives of / at v. 
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Lemma 3.6 For a smooth harmonic function f on M? whose second partial derivatives 
don't all vanish at any point, there is a bounded function f e on "Kt such that 

f e (z) = f(z) + e 2 f 2 (z) 

has discrete Laplacian of order 0(e 3 ), where f 2 (z) depends only the second derivatives 
of f at z and on the local structure of the graph "Kt oi z. 

Proof: We have exact formulas for one discrete harmonic function, the conjugate 
Green's function on elKy, and we know its asymptotics (Lemma 13. 5p . equation (JTT|) . 
which are 

G*(w, z 2 ) R* r)(zi,z 2 ) = —lm(clog(z 2 - z{j) 

for a constant c depending on the local structure of the graph near w, and where z\ is 
a point in 0(w). We'll let w be the origin in elK-r and z\ = 0; 77(0, z 2 ) is a continuous 
harmonic function of z 2 . 

For zfl 2 consider the second derivatives of the function /, which by hypothesis 
are not all zero. There is a point z 2 = (3(z) G M? at which n(0,z 2 ) has the same 
second derivatives. Indeed, / has three second partial derivatives, f X x,fxy, and f yy , 
but because / is harmonic f yy = —f xx - We have ^-^(0, z 2 ) = ^Im^-, and the second 
partial derivatives of r\ are the real and imaginary parts of const/ z 2 , which is surjective, 
in fact 2 to 1, as a mapping of M. 2 — {0} to itself. In particular there are two choices 
of z 2 for which f(z) and n{(5{z)) have the same second derivatives. Since / and n are 
smooth, by taking a consistent choice, (3 can be chosen to be a smooth function as well. 

Consider the function 

f e (z) = f(z) + GX(0,/3(z))- V (0,p(z)), 

where G* x is the discrete conjugate Green's function and A is chosen so that %r(A) (see 
section l3TT|) has local structure at (5{z) identical to that of "Kt at z. 

We claim that the discrete Laplacian of f e is 0(e 3 ). This is because G* x (0, (3(z)) is 
discrete harmonic, and f(z) — n(0,(3(z)) has vanishing second derivatives. 

We also have that G*(0,/3(z)) - r)(O,0(z)) = 0(e 2 ), see Lemma [33] above. □ 



3.4 Kn^ in constant-slope case 

Theorem 3.7 In the case of constant slope v and a bounded domain U, let £ be a 
conformal diffeomorphism from 4>(U) to H. When b and w are converging to different 
points as e — > we have 

K-i(h wl = 1 Tm ( g / (^(b))F(w)F(b) mbfg(b)\ 

^ ' > 2vrRe(F(w))F(b) \£(<j>(b)) - £fa(w)) + ^(0( b )) - ^R) / 

(18) 
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Proof: The function Gg T (w,v) is equal to the function (|16p for the whole plane, 
plus a harmonic function on St whose boundary values are the negative of the values 
of (fl6j) on the boundary of St- 

Since discrete harmonic functions on St are close to continuous harmonic functions 
on U, we can work with the corresponding continuous functions. 

From (|17p we have 

G5fr(w,t;) = ^Im ( Re ^ )} log(v - zx) )+es(z 1 ,v)+0(er, (19) 



2vr 



where zi is a point in face </>(w). The continuous harmonic function of v on {7 whose 
values on dU are the negative of the values of (fT9l) on <9C7 is 



1 t ( H w ) , / A 
^ Im (Re^R l0g(w -^J + 



(20) 

where S2 is smooth. 

The discrete Green's function Gg T must be the sum of (fT9j) and ([20j) : 

G^w.r,) = i-Im log(^) + i^log(^) +6.3(^1, t,)+0(e) 

Differentiating gives the result (as in Lemma I3.5H . □ 

It is instructive to compare the discrete conjugate Green's function in the above 
proof with the continuous conjugate Green's function on U which is 

1 1 



g*( Zl ,z 2 ) = — Imlog(£(zi) - t(z 2 )) + — Imlog(£(zi) - Z{z 2 )). 
3.4.1 Values near the diagonal 

Note that when b is within O(e) of w, and neither is close to the boundary, the discrete 
Green's function G*(w,v) for v on b is equal to the discrete Green's function on the 
plane G^ Kt (w,v) plus an error which is O(l) coming from the corrective term due to 
the boundary. The error is smooth plus oscillations of order 0(e 2 ), so that within 0(e) 
of w the error is a linear function plus 0(e). Therefore when we take derivatives 

^(b,w) = K^(b,w) + 0(l), 

which, since is of order 0(e _1 ) when |b — w| = 0(e), implies that the local 

statistics are given by \x v . 

Theorem 3.8 In the case of constant slope v , the local statistics at any point in the 
interior of U are given in the limit e — ► by \x v , the ergodic Gibbs measure on tilings 
of the plane of slope v. 
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4 General boundary conditions 



Here we consider the general setting: U is a smooth Jordan domain and h, the nor- 
malized asymptotic height function on U, is not necessarily linear. 

4.1 The complex height function 

The equation (pQ) implies that the form 

log(<J> - l)dx - log(^- - l)dy 

is closed. Since U is simply connected it is dH for a function H : U —* C which we call 
the complex height function. 

The imaginary part of H is related to h: we have arg(<3> — 1) = n — 9 C (Figure [2]) 
and arg(^j — 1) = 9 a — ir, which gives ImdH = (it — 9 c )dx + (ir — 9 a )dy. From ([3]) we 
have dh = (3p c — l)dx + (3p a — l)dy, so 

3 

-ImdH = 2{dx + dy) - dh. 

7T 

The real part of H is the logarithm of a special gauge function which we describe below. 
We have 

H & = log($-l) (21) 
H y = -log(i-l) (22) 

H ™ ~ $37 -^TT (23) 

H » = -f^lT- (24) 



4.2 Gauge transformation 

The mapping $ is a real analytic mapping from U to the upper half plane. It is an 
open mapping since Im(<]? y ) = — Im<l> ^ 0, but may have isolated critical points. 
The Ahlfors-Bers theorem gives us a diffeomorphism <p from U onto the upper half 
plane satisfying the Beltrami equation 

# d®. $ _ e ivr/3 
dz _ dz_ _ _ c 

# _ _ $ _ e -in/3 ' 
dz dz 

that is 0^ = Such a exists by the Ahlfors-Bers theorem [JJ. It follows that 

$ is of the form f(cj)) for some holomorphic function / from H into HI. Since dU is 
smooth, 4> is smooth up to and including the boundary, and <f>£, 4> y are both nonzero. 
For white vertices of 9 define 

F(w) = e^W ^( w )(l + eM(w)), (25) 
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where M (w) is any function which satisfies 
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The existence of such an M follows from the fact that the ratio of the coefficients 
of Mi and M y is -e H * +H * = $, so ([26]) is of the form 

Mi + = J(x, y) 

for some smooth function J. This is the 3 equation in coordinate 0. We don't need to 
know M explicitly; the final result is independent of M. We just need its existence to 
get better estimates on the error terms in Lemma [Jj] below. 
For black vertices b define 



F(b) = e-Wy07M(l-^H), (27) 
where w is the vertex adjacent to and left of b and M(w) is as above. 
Lemma 4.1 For each black vertex b with three neighbors in S we have 

^F(w)K(w,b) =0(e 3 ) 

w 

and for each white vertex w we have 

^K(w,b)F(b) = 0(e 3 ). 

b 

Proof: This is a calculation. Let w,w — ex,w + ey be the three neighbors of b. 
Then, setting H = H(w), (f) y = 0^(w), and M = M(w) we have 



F(w-ex) = e \{H-eH i+ ^H m -^H^O{^))^^ _ e ^ £ + + (e*)(l+eM-e 2 M & +0(e 3 )) 



F(w+ey) = e \{H+eH^H^i H m +0{*)) J ^ + ^ + + 0(e 3)( 1+eM+e 2 M??+0(e 3 ))> 

The sum of the leading order terms in F(w) + F(w — ex) + -F(ui + ey) is 



; ^70~(1 + e~ H * + e^) = e^70(l + + _ 



0. 



The sum of the terms of order e is fe^ \f4>y times 



J y yy 
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and the sum of the order-e 2 terms is e 2 e « 




4>y times 



M & + e^M# + e 



— #4 




+ 



) 



0. (28) 



A similar calculation holds at a white vertex, and we get the same expression for the 
e 2 contribution (changing the signs of M, H, d/dx and d/dy gives the same expression). 

□ 

It is clear from this proof that the error in the statement can be improved to any 
order 0(e 3+k ) by replacing M(w) with Mq(w) + eMi(w) + • • • + e k Mj t (\N), where each 
Mj satisfies an equation of the form (I26|) except with a different right hand side — the 
right-hand side will depend on derivatives of H, <f> and the Mj for i < j. For our proof 
below we need an error 0(e 4 ) and therefore the Mq and M\ terms, even though the 
final result will depend on neither Mq nor M\. 

4.3 Embedding 

Define a 1-form 



on edges of 9, where F is defined in (|25|27p (since K(w, b) = e is a constant, this is an 
unimportant factor for now, but in a moment we will perturb K). By the comments 
after the proof of Lemma l4.lt the dual form Q* on 9* can be chosen to be closed up to 
0(e 4 ) and so there is a function cj) on 9*, defined up to an additive constant, satisfying 
# = fl* + 0(e 3 ). 

In fact up to the choice of the additive constant, <p is equal to <fi plus an oscillating 
function. This can be seen as follows. For a horizontal edge wb we have F(w)F(b) = 
4>y(w) + 0(e 2 ). Thus on a vertical column {wibi, . . . , w^b^} of horizontal edges we have 

k k k 



Y, tt*(wA) = e ^CF(wi) + F(yn))F(bi) = e ^ <^( Wi ) + e ^ F( Wi )F(bi) + 0(e 3 ). 



The first sum gives the change in <f> from one endpoint of the column to the other, 
and the second sum is oscillating (F(w{) and F(b{) have the same argument which is 
in (0, 7r) and which is a continuous function of the position) and so contributes 0(e). 
Similarly, in the other two lattice directions the sum is given by the change in plus 
an oscillating term. 

Therefore by choosing the additive constant appropriately, (j> maps 9* to a small 
neighborhood of H in the spherical metric on C, which shrinks to EI as e — > 0. 



n(wb) = 2Re(F(w)) J F(b)ET(w, b) 



i=i 



i=l 



i=l 
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The map <j) has the following additional properties. The image of the three edges 
of a black face of 9* are nearly collinear, and the image of a white triangle is a triangle 
nearly similar to the a, b, c-triangle, and of the same orientation. Thus it is nearly a 
mapping onto a T-graph. In fact near a point where the relative weights are a, b, c 
(weights which are slowly varying on the scale of the lattice) the map is up to small 
errors the map of section [3j 

We can adjust the mapping (f) by 0(e 3 ) so that the image of each black triangle 
is an exact line segment: this can be arranged by choosing for each black face a line 
such that the (/>-image of the corresponding black face is within 0(e 3 ) of that line; the 
intersections of these lines can then be used to define a new mapping ^ : 9* —> C which 
is an exact T-graph mapping. 

The mapping ^ will then correspond to the above 1-form Q but for a matrix K 
with slightly different edge weights. Let us check how much the weights differ from 
the original weights (which are e). As long as the triangular faces are of size of order e 
(which they are typically), the adjustment will change the edge lengths locally by 0(e 3 ) 
and therefore their relative lengths by l + 0(e 2 ). There will be some isolated triangular 
faces, however, which will be smaller — of order 0(e 2 ) because of the possibility that 
T(w) might be nearly pure imaginary. We can deal with these as follows. Once we have 
readjusted the "large" triangular faces we have an exact T-graph mapping on most of 
the graph. We can then multiply T(w) by A = i and F(b) by A = —i: the readjusted 
weights give (for most of the graph) a new exact T-graph mapping (because we now 
have KF = FK = exactly for these weights), but now all faces which were too small 
before become O(e) in size and we can readjust their dimensions locally by a factor 
l + 0(e 2 ). 

In the end we have an exact mapping ^ of 9* onto a T-graph 9t and it distorts 
the edge weights of K by at most 1 + 0(e 2 ). We shall see in section [5] that this is close 
enough to get a good approximation to K^ 1 . 

In conclusion the Kasteleyn matrix for 9d is equal to 2Re(T(w))T(b)i^(w, b) where 
K has edge weights e + 0(e 3 ). 

4.4 Boundary 

Along the boundary we claim that the normalized height function of 9d follows u. Since 
9 d arises from a T-graph, we can use its canonical flow (section I3.2.2P . Near any given 
point the canonical flow looks like the canonical flow in the const ant- weight case — since 
the weights vary continuously, they vary slowly at the scale e, the scale of the graph. 
Since the canonical flow defines the slope of the normalized height function, we have 
pointwise convergence of the derivative of h along the boundary to the derivative of u. 
Thus h converges to u. In fact this argument shows that the normalized height function 
of the canonical flow converges to the asymptotic height function in the interior of U 
as well. 
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5 Continuity of K 1 

In this section we show how K~ l changes under a small change in edge weights. 

Lemma 5.1 Suppose < 5 -C e. If 5' is a graph identical to Sd but with edge weights 
which differ by a factor 1 + 0(5), then 

K^=K^(l + 0(8/e)). 

In particular since 5 = e 2 in our case this will be sufficient to approximate K^ 1 to 
within 1 + 0(e). 

Proof: For any matrix A we have 

(oo 
i + ^(-iy(8ey(AK- 1 y 

as long as this sum converges. 

From Theorem 16.11 below we have that iCg, 1 (b,w) = 0(l/|b — w|). When A rep- 
resents a bounded, weighted adjacency matrix of S_d, the matrix norm of AK^ 1 is 
then 

\\AK- l \\ < maxV I V A(W ,b)K~ 1 (b,\N)\ < 3amax V — - — r = 0(e~ 2 ) 

w b w^w' 

where a is the maximum entry of A. In particular when 8 <C e we have 

(K + SeA)' 1 = K~ l (l + 0(8/ e)). 

□ 



6 Asymptotic coupling function 

We define Kg D as in (fT3|) using the values (|25l27p for F (and Kg is the adjacency 
matrix of S). 

Theorem 6.1 // b, w are not within o(l) of the boundary ofU, we have 



When one or both o/b,w are near £/ie boundary, but they are not within o(l) o/ eac/i 
of/ier, tfien Kg^(b,w) = 0(1). 

Proof: The proof is identical to the proof in the constant-slope case, see Theorem 
13.71 except that £(f> there is <f> here. The £' factors from (I18p are here absorbed in the 
definition of the functions <p and F. □ 

As in the case of constant slope, when b and w are close to each other (within 
0(e)) and not within o(l) of the boundary, Theorem 13. 81 applies to show that the local 
statistics are give by fi u . 
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7 Free field moments 



As before let cf> be a diffeomorphism from U to HI satisfying <f>± = —&(f)y where $ is 
denned from h as in section 14.11 

Theorem 7.1 Let h be the asymptotic height function on Sd, h e the normalized height 
function of a random tiling, and h = ^^-(h e — h). Then h converges weakly as e — > 
to the pull-back under (ft of "5 , the Gaussian free field on H. 

Here weak convergence means that for any smooth test function tp on U, zero on 
the boundary, we have 



e 2 £>(/)&(/) 

So 



^(x)3{(i)(x))\dx\ 2 , 



where the sum on the left is over faces / of Sd- 

Proof: We compute the moments of 3~. Let tpi, . . . , ipi- be smooth functions on U, 
each zero on the boundary. We have 



E 



E Mfi)Hfi) ) ■ ■ U 2 Mfk)Hfh 



= e 2k £ Mfi)---Mfk)nh(fi)...h(f k )}. 

fl,—tfk 

From Theorem 17.21 below the sum becomes 

ETO(xi)) . . . H(t>(x k ))] Yl ^(x^dxtf + o(l), 



J7 JU 



that is, the moments of h converge to the moments of the free field <j)*(9'). Since the 
free field is a Gaussian process, it is determined by its moments. This completes the 
proof. □ 

Theorem 7.2 Let s\, . . . , s& G U be distinct points in the interior of U. For each e 
let fi, $2, ■ ■ ■ , fk be faces of ' 9d, with fi converging to Si as e — > 0. If ' k is odd we have 

limE[h(f 1 )...h(f k )]=0 

e— >0 



and if k is even we have 



k/2 



KmE[h(fi)...h(f k )}= n G ^(Mjy-i))'^(M2i))) 

pairings a j=l 



where 



G(z,z') 



2tt 



log 



z — z 



z — z' 
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is the Dirichlet Green's function on HI and the sum is over all pairings of the indices. 
If two or more of the Si are equal, we have 

E[h(f 1 )...h(f k )]=0(e- £ ), 

where i is the number of coincidences (i.e. k — £ is the number of distinct Si). 

Proof: We first deal with the case that the Sj are distinct. Let 71, . . . , 7^ be pairwise 
disjoint paths of faces from points s[ on the boundary to the /j. We assume that these 
paths are far apart from each other (that is, as e — * they converge to disjoint paths). 
The height h e (fi) can be measured as a sum along 7$. 

We suppose without loss of generality that each 7$ is a polygonal path consisting 
of a bounded number of straight segments which are parallel to the lattice directions 
x,y,z. In this case, by additivity of the height change along 7$ and linearity of the 
moment in each index, we may as well assume that 7$ is in a single lattice direction. 

Now the change in h along 7, is given by the sum of a^- — E(a^) where a^- is the 
indicator function of the jth edge crossing 7$ (with a sign according the the direction 
of 7i). So the moment is 

nkh) ■ ■ ■ kfk)] = E E [Ki - E M) • • • k* - Efoyj)]. 

Jl671i— ,3k£lk 

If Wij i} bij i are the vertices of edge a^, this moment becomes (see [8]) 

if- 1 ^,™^) ... K-^bfc^.wi^) 

K _1 ( b lJl 1 w kj k ) ■■■ K ~ 1 ( b k-lj k _ 1 . w fcj fc ) 



n n 



In particular, the effect of subtracting off the mean values of the is equivalent to 
cancelling the diagonal terms K~ 1 (bij i ,Wij i ) in the matrix. 

Expand the determinant as a sum over the symmetric group. For a given per- 
mutation <7, which must be fixed-point free or else the term is zero, we expand out 
the corresponding product, and sum along the paths. For example if a is the /c-cycle 
a = (12 . . . k), the corresponding term is 

sgn(ff) E (n^^' 6 ^)) ^H b iii' w 2i2)^ _1 ( b 2j 2 ,w3i 3 )---i ; s:~ 1 ( b fci fe ,wi jl ) = 

jl,-,jh \i=l I 



_ -l — ) k V ( F(bl ^ )F( ' N2 ^ ) F(bl 31 )F(™ 2]2 ) F(bi 31 )F(w 2j2 ) F(b ln )F(w 2j2 ) 



Jl.— .Jfc 



<l>(blj 1 ) - 4>{™2j 2 ) ^(bijj) - 0(w 2 j 2 ) ^"(blj-j) - ^(w 2 j 2 ) <H b lii) - 0(w 2 j 2 ) 



f (bfcj fc )J'('"lj 1 ) f(b fcjfc )-F(w ljl ) F(b fcjfc )F(w 1;| - 1 ) f (bfcj fc )f(wij 1 ) 



■XblJ! ) - ^( w 2j 2 ) 0C b fcJ fc ) - ^( w lJl) 0C b fcJ fc ) - ^( w lJl) 0C b fcJ fc ) - ^( w lJl) 



(30) 



plus lower-order terms. Multiplying out this product, all 4 fc terms have an oscillating 
coefficient (as some ji varies) except for the terms in which the pairs F(bjjJ and F(vjij i ) 
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in the numerator are either both conjugated or both unconjugated for each i. That 
is, any term involving F{bij i )F(\Nij i ) or its conjugate will oscillate as ji varies and so 
contribute negligibly to the sum. There are only 2 k terms which survive. 
Let Zi denote the point z% = 4>(bij z ) ~ <^(wy\). 

For a term with F(byJ and F(wij i ) both conjugated, the coefficient 



is equal to dzj, otherwise it is equal to dzi, where dzi is the amount that 4> changes 
when moving by one step along path 7$, that is, when ji increases by 1. 
So the above term for the fe-cycle a becomes 



e n 




dz[ £l) . . . dz^ k) 



' 4 ™ £1 ,..ir= ± i Vt ) J ^ ' ' ' J ^ (4 £l) - 4 £2) )(4 £2) - 4 £3) ) • • • (4 £fc) - 4 £l) ) 



plus an error of lower order, where zj 1 ^ = Zj and Zj ^ = Zj, with similar expressions 
for other a. 

When we now sum over all permutations a, only the fixed-point free involutions do 
not cancel: 

Lemma 7.3 (|2j) For n > 2 let C n be the set of n- cycles in the symmetric group S n . 
Then 

n 1 
where the indices are taken cyclically. 

Proof: This is true for n odd by antisymmetry (pair each cycle with its inverse). 
For n even, the left-hand side is a symmetric rational function whose denominator 
is the Vandermonde n^-Ozj — z j) aim whose numerator is of lower degree than the 
denominator. Since the denominator is antisymmetric, the numerator must be as well. 
But the only antisymmetric polynomial of lower degree than the Vandermonde is 0. 

□ 

By the lemma, in the big determinant all terms cancel except those for which 
a is a fixed-point free involution. It remains to evaluate what happens for a single 
transposition, since a general fixed-point free involution a will be a disjoint product of 
these: 

/■<H s l) f<P(s2) dz\dz*2 /'<M s l) f4 l (s2) dz\dz2 f&(. s l) /'<M S 2) dz\dz~2 f$( s l) f<P(s2) dz\dz 2 



(4tt 



rf('2) dz 1 dz 2 f»('l) r<t>(s2) dzidz2 f*( = l) / , «(s 2 ) dz 1 dz 2 r<t>(si) r<t>(s2> dz 

1 dz x dz 2 



( 4 ™) 2 J^j (zi - z 2 ) 

= i-G(0( Sl ),^ 2 )) 



2 



where we used d>(s' 4 ) £ 
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Now suppose that some of the Sj coincide. We choose paths 7$ as before but suppose 
that the j{ are close only at those endpoints where the s, coincide. Let 5 = y/e. For 
pairs of edges 0*^,0^., on different paths, both within 5 of such an endpoint we use 
the bound JC _1 (b,w) = O(^), so that the big determinant, multiplied by the prefactor 
Y\ { K{\Nij^ bijj, is 0(1). In the sum over paths the net contribution for each coincidence 
is then 0(5/e) 2 = 0(e~ 1 ), since this is the number of terms in which both edges of two 
different paths are near the endpoint. □ 



8 Boxed plane partition example 

For the boxed plane partition, whose hexagon has vertices 

{(1,0), (1, 1), (0, 1), (-1, 0), (-1, -1), (0, -1)} 

in (x,y) coordinates, see Figure O it is shown in [11] that satisfies (— <3?x + y) 2 = 
1 - + $ 2 , or 

1 - 2xy + v^x 2 - xy + y 2 ) - 3 

= 2(T=P) 

for (x,y) inside the circle x 2 — xy + y 2 < 3/4. $ maps the region inside the inscribed 
circle with degree 2 onto the upper half-plane, with critical point (0, 0) mapping to 
e 7 ™/ 3 . If we map the half plane to the unit disk with the mapping z i— ► z ~ e ^.J 3 , the 
composition is 

(31) 



$(x,y) - e"/ 3 2r 2 - 3 + V9 - 12r" 



$(x,y) - e~ i7r / 3 2(y-tox) 2 

(where r 2 = x 2 — xy + y 2 ) which maps circles concentric about the origin to circles 
concentric about the origin. To see this, note that z = i(y — lox) defines the standard 
conformal structure, and \z < \ 2 = x 2 — xy + y 2 = r 2 , so that the right-hand side of the 
equation (j3Tj) is f{r)/z 2 . 

Note also that the Beltrami differential of <£, which is 



<J> _ g-TTi/3 ' 

satisfies 

A» = — > 

that is, it is its own Beltrami differential! This is simply a restatement of the PDE ([1]) 
in terms of 

The diffeomorphism (f> from the region inside the inscribed circle to the unit disk is 
also very simple, it is just <fi = yfji, or 

2r 
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The inverse of 4> which maps D to U is even simpler: it is 



<T\z) 



l -h U| 



This map can be viewed as the orthogonal projection of a hemisphere onto the plane 
through its equator, if we identify D conformally with the upper hemisphere sending 
to the north pole. 

In conclusion if the domain U is the disk {(x,y) | x 2 — xy + y 2 < r 2 }, where 
r 2 < 3/4, and the height function on the boundary of U is given by the height function 
h of the BPP on dU, then the h on U will equal the h on BPP restricted to U, and 
the fluctuations of the height function are the pull-back of the Gaussian free field on 
the disk of radius ^~H~ Arl under <j>. 
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